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ABSTRACT 

Using two archival XMM-Newton observations, we identify strong X-ray flux variations in NGC 
6946 X-l indicating it is the most variable ultraluminous X-ray source (ULX) on mHz time scales 
known so far. The 1-10 keV lightcurve exhibits variability with a fractional rms amplitude of 60% 
integrated in the frequency range of 1-100 mHz. The power spectral density of the source shows a 
flat-topped spectrum that breaks at about 3 mHz with possible quasi-periodic oscillations (QPOs) 
near 8.5 mHz. Black hole binaries usually produce strong fast variability in the hard or intermediate 
state. The energy spectrum of NGC 6946 X-l is dominated by two components, a 0.18 keV thermal 
disk and a power law with a photon index of «2.2, which is consistent with the intermediate state. 
The characteristic time scales of the X-ray emission suggests that the ULX may contain a black hole 
with a mass on the order of 10 3 solar masses. 

Subject headings: black hole physics — accretion, accretion disks — X-rays: binaries — X-rays: 
individual (NGC 6946 X-l) 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are bright, non- 
nuclear X-ray sources found in nearby galaxies with lu- 
minosities above 3 x 10 erg s _1 assuming isotropic 
emission, which is higher than the Eddington limit of 
a 20 Mq black hole. Their nature remains a ma- 
jor astrophysical puzzle. They may be intermediate 
mass black holes of 10 2 - 10 4 M Q ([Colbert fc Mushotzkvl 
1999), or dinary stellar mass black holes of ~10 M Q with 
beamed (jKing et al.ll200lHKording et al.ll2002h or super- 
Eddington radiation (|Begelmanll2006f) . or, massive stellar 
black holes of <100 Mq formed in low metallicity envi- 
ronm ent with a high accretion rate (jZampieri fc Roberts! 
2009). It is also possible that the high apparent lumi- 
nosity i s due to a combination of two or three of these 
factors (jPoutanen et alJ 120071 : King! 120081 ). Due to the 
large distance and the brightness of the disk, dynamical 
measurement of the compact object mass via the motion 
of the companion is difficult and has not yet been accom- 
plished for any ULX. X-ray spectroscopy could shed light 
on the black hole mass, but is only reliable in the thermal 
dominant state, which is rarely found in ULXs excep t 
one case firmly identified in M82 (|Feng fc Kaaretll2010D . 
Another means to constrain compact object masses is X- 
ray timing. Characteristic time scales of X-ray emission, 
such as quasi-periodic oscillations (QPOs) and the break 
frequency in broad-band timing noise, have proved useful 
in constraining the mass of compact objects via a model 
independent calibration (cf. McHardy ct al. 2006). 

Low-frequency QPO s have been detected in thre e 
ULXs, M82 X41.4 +60 dStrohmaver fc Mushotzky| [2003T), 
NGC 5408 X-l (IStrohmayer et all 120071). and M82 
X42.3+59 (|Feng et all l2010f ). To begin to understand 
these QPOs, an analogy with the QPOs detected from 
stellar mass black holes may be helpful. Low-frequency 
QPOs are classified into three main classes in stellar mass 
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black holes, as type A, B, or C based on their amplitude, 
coherence, ha rmonics, red noise, etc. (for a thorough 
discussion see iCasella et all [20051 ) . The QPOs in M82 
X41.4+60 and NGC 5408 X-l are strong and narrow, 
and ap pear above strong flat-topped, band-limited noise 
(BLN) dStrohmaver fc M ushotzkvl 120031 : iFeng fc Kaaretl 
120071: IStrohmayer et all 120071) . The QPO frequencies 
vary in a wide range of 50-100 mHz for M82 X41.4+60 
and 10-20 mHz for NGC 5408 X-l and the BLN shows a 
break from a power-law index of to —1 near the QPO 
frequency (iDewangan et ap2006l : lMucciarelli et al.ll2006t 
IStrohmayer fc M ushotzkv 2009). These properties are 
typical of the type C QPOs, which are the most com- 
mon low-frequency QPOs found in Galactic black hole 
binaries usually appearing in the hard and intermedi- 
ated states, and occasionally in t he steep power-law state 
(|McClintock fc Remillardl 120061 ) . Recently, QPOs were 
discovered in another ULX in M82, X42.3+59. These 
QPOs show different characteristics: they are strong 
and broad, appear with very weak or absent red noise, 
detected only when the source is more luminous than 
10 erg s , and their central frequency varies in a nar- 
row range from 3 to 4 mHz (jFeng et all 120101 ). These 
properties are unlike from type-C QPOs, but instead, 
are similar to those of type A/B QPOs, which show 
up in the steep power-law state. However, a survey of 
ULXs impli es that the majority do not show strong fast 
variability (IFeng fc Kaaretl 12 005). From a sample of 16 
bright ULXs, iHeil et aLl (|2009l ) concluded that in some 
ULXs the variability is significantly suppressed compared 
to bright black hole binaries and active galactic nuclei. 
Thus, detection of characteristic time scales in ULXs 
provides clues to their masses, but appears only rarely. 
Therefore, it is of significant interest to search for fast 
variability from ULXs. Here, we report the detection of 
strong X-ray variability in the mHz range in NGC 6946 
X-l. 

NGC 6946 is a n early face-on spiral galax y at a dis- 
tance of 5.9 Mpc (jKarachentsev et al.l [2000h . The X- 
ray s ource NGC 6946 X-l was fir st detected with Ein- 
stein ([Fabbian o fc Trinc hieri 1983), and later was found 
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to be spatially associated with a nebula commonly re- 
ferred to_as_^lF_JJ_j ; ndJhojaglit to be a supernova rem- 
nant dMatonkk k Fesed 119971) on the basis of R OSAT 
observations (ISchlegell 119941 IBlair k Feseri Il99l . Us- 
ing Chandra observations. IRoberts k Colbertl ( |2003l ) de- 
tected X-ray variability on time scales of minutes and 
argued that the X-rays arise from a black hole X-ray 
binary inside the nebula rather than extended emission 
from the nebula itself. The compact nature of the X- 
ray source was further strengthened with the detection 
of ran dom variations in the X-ray flux on long time 
scales (|Fridriksson et al.l 120081 ) . The surrounding neb- 
ula MF 16 is thought to be a product of shock ioniza- 
tion ([Blair et al.l 120011) plus unusually strong photoion- 
izatio n possibly originating fr o m the central c o mpac t ob- 
ject (|Abolmasov et alJl2008| ). Kaaret et all (|201Q[ ) de- 
tected an ultraluminous UV counterpart to the X-ray 
source right in the center of MF 16, which confirmed 
this scenario. Here, using two XMM-Newton observa- 
tions available in the archive, we analyzed its timing and 
spectral behavior in § [2j and discussed the constraint on 
its physical nature in § [3] 

2. OBSERVATIONS AND DATA ANALYSIS 

We examined all archival XMM-Newton observations 
of NGC 6946 with exposures longer than 20 ks to search 
for timing noise of X-l at low frequencies. Among 
these observations, four out of six are not suited for this 
work because of high particle background, with ~30% 
background contamination in the source region. The 
remaining two observations (ObsIDs 0500730201 and 
0500730101) were adopted for spectral and timing anal- 
ysis. 

In these two observations, both PN and MOS were 
operated in full frame mode. We used SAS 8.0.1 and cal- 
ibration files current as of April 2010 for standard data 
reduction. For both the PN and MOS cameras, each 
energy spectrum was extracted from a 12"-radius circu- 
lar region to avoid the chip gap, and background from a 
nearby, circular, source-free region was subtracted. We 
fitted the spectra with th e best model as discussed in 
IRoberts k Colbert] (|2003l ): a cool multicolor disk plus 
a power-law component. All spectral parameters with 
errors are shown in Table [TJ The power-law component 
contributes ^50% to the total flux in the 0.3-10 keV band 
and ^80% in the 1-10 keV band. The two observations 
were made within a week and their spectral parameters 
are consistent within the 90% errors. In particular, if 
we impose the same absorption column density on the 
two fits, almost identical results are obtained. Thus, we 
claim no significant spectral variation between the two 
observations, and, hence, we simultaneously fitted the 
two observations with a single model, as above, with the 
addition of a constant multiplying the flux normalization 
for ObsID 0500730101. As shown in the fourth column 
of Tabled! this model fits well, with x 2 /dof = 516.9/511. 
The constant is estimated to be 0.98 ± 0.03, consistent 
with 1 within the 90% error, which indicates an approx- 
imately constant flux. 

Due to the large effective area and fine time resolution 
of PN, we created a light curve using PN events in the 
GTIs with a time step equal to the frame time (~ 73.4 
ms) . The PN lightcurves shown in Figure Q] have a time 
step of 200 s for visual inspection only. A Fourier trans- 



Table 1 

Spectral Fits to the XMM-Newton Spectra 



ObsID 


0500730201 


0500730101 


Combined 


Observation date 


2007-11-02 


2007-11-08 




Exposure (ks) 


37 


33 




jV H (10 22 cm' 2 ) 


0.35 ±0.05 


0.29 ± 0.05 


0.32 ±0.03 


T in (keV) 


0.17 ±0.02 


0.20 ±0.03 


0.18 ±0.02 


R in (10 3 km) 


7 o+5.6 
—2 3 


4 2+ 3 ' 3 
^• z -1.3 


c c+2.7 
°-°-1.4 


r 


2.33 ±0.13 


2.05 ±0.17 


2.19 ±0.10 


fx (IO" 13 erg/cm 2 /s) 


8.3 ±0.3 


8.6 ±0.3 


8.5 ±0.2 


Lx (10 39 erg/s) 


11 7+ 13 

1±.<_ 2 7 


o 4+3.3 


10.0 ± 1.6 


L d (10 39 erg/s) 
X 2 /dof 


5 9+ 3 ' 2 


42 +i.8 


5.0 ± 1.4 


289.9/290 


221.5/217 


516.9/511 



Note. — The combined spectrum is described in the text. Nu 
is the absorption column density, T; n is the temperature of inner 
disk, R i n is the radius of the inner disk assuming a zero inclination 
angle, T is the power-law photon index, fx is the observed flux 
in 0.3-10 keV, Lx ls the luminosity corrected for absorption in 
0.3-10 keV, L d is the disk luminosity corrected for absorption in 
0.3-10 keV. Errors arc quoted at a confidence of 90%. 




Figure 1. X-ray lightcurve of NGC 6946 X-l from PN data in the 
energy range of 1-10 keV for the observation 0500730201 (top) and 
0500730101 (bottom). There are three GTIs in 0500730201 and two 
in 0500730101. The PSDs are created from the GTIs only. 



form was performed within each GTI with 2 15 points 
and individual spectra were averaged to create the final 
one. We grouped the frequency channels linearly by a 
factor of 2 followed by a logarithmic binning by a factor 
of 10 1 / 30 . These two steps led to at least 26 and 20 power 
points for observations 0500730201 and 0500730101, re- 
spectively, combined into a single frequency bin as re- 
quired to model the powe r spectrum using minimum y 2 
fitting flvan der Klisl 1198 9: Pap adakis k Lawrence! 119931 : 
iHeil et al.ll2009f >. The PSDs for observations 0500730201 
and 0500730101 from photons in the energy range of 0.3- 
10 keV are shown in Fig ure and Eb, resp ectively, and 
are normalized following iLeahv et al.l (|1983l ) . They both 
show strong variability at low frequencies above the Pois- 
son level, which is a constant value of 2 in the Leahy 
normalization. As for the energy spectra, we combined 
the PN power spectra from the two observations, see Fig- 
ure [5J;. 

To quantify the variability we fitted the power spectra 
using XSPEC vl2. First we fitted the 0.3-10 keV PSD 
with a power-law model for the red noise plus a con- 



Low Frequency Variability in NGC 6946 X-l 



3 




10 

Frequency (Hz) 

Figure 2. X-ray power spectra of NGC 6946 X-l calculated from 
the PN data, (a) PSD of observation 0500730201; (b) PSD of ob- 
servation 0500730101; (c) PSD of the two observations; (d) similar 
to (c) but in 1-10 keV; The PSDs of (a)-(c) are vertically shifted for 
clarity by 20, 12, and 6, respectively. The dashed lines indicate the 
fits with the two Lorentzian plus constant model, with parameters 
listed in Table [2] 

stant component for the Poisson noise, which gives a % 2 
of 142.5 with 88 degrees of freedom (dof). We then re- 
place the power-law component with a broken power- law 
with a fixed index of zero below the break (a flat-topped 
model). This improves the fit to x 2 /dof = 101.1/87, 
corresponding to a chance probability of 5 x 10~ 8 for 
the break as calculated using the F-test. A broad, 
zero-centered Lorentzian mod el is often used to char- 
acterize the flat-topped BLN ()Belloni et al.ll2002D . Fit- 
ting with a zero-centered Lorentzian plus a constant re- 
sults in a fit much better than the unbroken power-law, 
X 2 /dof = 98.8/88, again confirming the presence of a 



break. Adding a narrow Lorentzian to account for excess 
noise around 6-10 mhz, possible QPOs, further improves 
the fit to x 2 /dof = 77.7/85, suggesting a positive detec- 
tion with a chance probability ofl.3xl0~ 4 . There are 91 
independent frequency bins and the QPOs spreads over 
~ 5 bins, giving about 91/5 trials. Taking into account 
the number of trails, the chance probability is 2.3 x 10 -3 
equivalent to a significance of 3a. We note that this sig- 
nificance may be an underestimate, since the number of 
trails estimated in this way is best suited for periodic 
signals that concentrate all power in a single, indepen- 
dent frequency bin. Therefore, we conclude that we find 
a firm evidence for a break and a possible signature of 
QPOs in the power spectrum from NGC 6946 X-l. All 
timing parameters derived from the double Lorentzian 
model are listed in Tabled 

PSDs at three independent energy bands, 0.3-1 keV, 
1-2 keV, and 2-10 keV, were also produced following the 
same procedure, and the integrated rms in the frequency 
range of 1-100 mHz is 11%±3%, 53%±3%, and 70%±6%, 
respectively. This is similar to the behavior of black hole 
binaries where the timing noise amplitude often increases 
with energy. Since the variability appears mostly above 
1 keV, we produced a PSD in the 1-10 keV energy range, 
as shown in Figure However, the strength of the 
possible QPOs does not appear to increase with energy. 
New, long observations are needed to test the significance 
of the QPO. 

3. DISCUSSION 

Two XMM- Newton observations of NGC 6946 X-l re- 
veal strong X-ray variability in the energy range of 1- 
10 keV at short time scales with a band-limited spectrum 
that breaks at about 3 mHz and has a fractional rms am- 
plitude of 60% integrated over the range 1-100 mHz. We 
find evidence for possible QPOs with a central frequency 
of ~ 8.5 mHz. The detection of strong fast variability 
confirms the compact nature of the source, and, together 
with the high luminosity, suggests that the ULX contains 
an accreting black hole. 

A timing survey of bright ULXs indicates that only 
a small number of ULXs exhi bit strong variab ility at 
time scales shorter than 10 4 s pen et al.ll2009fl . NGC 
5408 X-l and M82 X41.4+60 are the typical two. 
The integrated rms amplitude in 1-100 mHz range in 
the 1-8 keV band is around 50 % for NGC 5408 X-l 
(|Strohmaver fc Mushotzkvl 120091) and slightly lower for 
M82 X41.4+60 (|Feng fc Kaaretj |2007l ). The integrated 
rms of NGC 6946 X-l in the same energy and frequency 
range is 59% ± 2%. Thus, NGC 6946 X-l becomes the 
most noisy ULX reported so far. 

The timing noise amplitude of black hole binaries 
varies dramatically in response to the emission state. 
When thermal emission from the accretion disk dom- 
inates the energy spectrum, timing noise is largely 
suppressed. In contrast, strong variability is always 
expected in the hard state, when the Comptoniza- 
tion component is dominant, and often appears, with 
a relatively lower amplitude, in the steep power-law 
state when the Comptonization co mponent is impor- 
tant (|McClintock fc Remillardl [2006) . This is also seen 
in ULXs. When M82 X41.4+60 changed its spectrum 
from a featureless power-law to a curved shape consis- 
tent with emission from a thermal disk, the previously 
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Table 2 

Timing parameters of NGC 6946 X-l fitted with a double Lorentzian model 





PSD in 0.3-10 kcV 




PSD in 1-10 kcV 




Flat-topped Noise 


QPO 


Flat-topped Noise 


QPO 


v (mHz) 





8.5 ±0.5 





8.3 ±0.4 


FWHM (mHz) 


4.9 ± 1.1 


3.8 ±1.5 


6.8 ± 1.6 


3.1 ± 1.6 


Norm (10 -2 ) 


4.1 ±0.7 


1.1 ±0.3 


5.3 ±0.9 


1.0 ±0.4 


"max (mHz) 


2.5 ± 0.6 


8.7 ±0.7 


3.4 ±0.8 


8.4 ±0.5 


rms/mean (%) 


30.4 ±0.9 


26.4 ±3.9 


49.6 ± 1.4 


34.4 ±7.2 



Possion level 1.988 ± 0.003 1.993 ± 0.003 

R T (counts s" 1 ) 0.1712 0.0907 

R B (counts s" 1 ) 0.0033 0.0021 

Note. — u, FWHM, and Norm are the centroid, width, and integrated area of a Lorentzian model in XSP EC, respectiv ely. The units 
of 'Norm' are Leahy power multiplied by Hz. v max is the characteristic frequency where vP v peaks (Belloni ct al. 2003). The rms of 
flat-topped noise is calculated in 1-100 mHz frequency range. i?x and Rb are the total and background count rate in the source extraction 
region, respectively. Errors are quoted at a confidence level of 68.3% (If). 



detec ted QPOs and powe r continuum disappeared as 
well (iFeng fc Kaaretl [20101) . M82 X41.4+60 is believed 
to be in the hard state when the timi ng noise is seen 
based on its hard power-law spectrum ([Feng fc Kaaretl 
HoEj ), while NGC 5408 X-l is suggested to be in the 
steep power-law state due to its soft spectrum and the 
spect ral-temporal correlation ([Strohmaver fc Mushotzkvl 

The extremely high variability of NGC 6946 X-l sug- 
gests that the corona/Comptonization emission is very 
important and should be the dominant component sim- 
ilar to the hard or steep power-law state. With a frac- 
tional rms of 60%, it is more likely to be in the hard 
or intermediate state rather than the steep power-law 
state. The energy spectrum of NGC 6946 X-l from the 
two XMM-Newton observations is consistent with this 
scenario: the power-law component is dominant with a 
contribution of about 80% in 1-10 kcV and the power-law 
photon index of 2.2 ± 0.1 is intermediate between those 
typical of the hard and steep power-law states. How- 
ever, ULXs may not follow the general evolution pat- 
tern that most Galactic sources have. For example, GX 
339-4 shows a correlation between the variation ampli- 
tude and the spectral hardness, and the variation am- 
plitude exceeds 0.3 only when the source is hardest dur- 
ing the 'right branc h' on its hardness-intensity diagram 
([Belloni et al.l l2005). It is obviously difficult to interpret 
the energy spectrum of NGC 6946 X-l as being extremely 
hard. Conversely, GRS 1915+105 has shown an irregular 
pattern between the hardness and the total rms ampli- 
tude. The source often exhibits strong vari a tion up to a 
fracti onal rms of 70% ([Morgan et al.l 119971 : iMuno et al.l 
|1999T >. In this sense, NGC 6946 X-l is perhaps a better 
analog to GRS 1915+105. 

The break frequency of the broad-band noise is thought 
to be invers ely proportional to the mass of the compact 
objects (e.g. iMarkowitz et al.ll2003l) .Cvgnus X-l shows a 
similar break in the broad-band noise, from an index of 
to —1, at a frequency that varies i n the range of 0.04- 
0.4 H z depending on emission states ([Belloni fc Hasingerl 
|1990T ). The break frequency for NGC 6946 X-l is about 
3 mHz, which is 13-130 times lower than that in Cygnus 
X-l. This suggests that NGC 6946 X-l harbors a 130- 
1300 Mq black hole, as suming a compact o bject mass 
of 10M Q in Cygnus X-l ([Herrero et al.lll995| ). This esti- 



mate depends on the assumption that the low frequency 
break in the PSD scales linearly with black hole mass 
and extends well above 10M Q . 

The break frequency of broad-band noise is found to be 
tightly c orrelated with the QPO frequen cy in black hole 
binaries (jWiinands fc van der Klislll999D . The break and 
QPO frequencies reported here for NGC 6946 X-l appear 
to lie on a low frequency extension of the correlation mea- 
sured for stellar mass black holes, suggesting that simi- 
lar mechanisms produce the low frequency variability on 
both Galactic black hole binaries and this ULX. On the 
basis of the high QPO amplitude and the presence of the 
flat- top ped noise componen t, we classify the QPOs as of 
type C (|Casella et al. 2005). An interesting correlation 
has been found between the frequency of type-C QPOs 
and the photon index of the power-law component in 
the energy spectrum, and it is suggested that the QPO 
frequency is inversely scale d with the black hole mass 
at a given photon index (Shaposh nikov fc Titarchuckl 
12009ft . For a photon index range of 2.1-2.3, the QPO 
freq uency of GRS 1915+ 105 varies roughly between 1-2 
Hz ([Vignarca et al.ll2003[ ). which is 120-240 times higher 
than the frequency of the possible QPOs in NGC 6946 
X-l. Thus, assuming a positive detection of type-C 
QPOs in N GC 6946 X-l and a 14 M p black hole in GRS 
1915+105 (|Harlaftis fc Greineri[200l . we derive a black 
hole mass of (1 - 4) x 10 3 M for NGC 6946 X-l. We 
caution that the QPO frequency-photon index relation 
has been verified only for stellar mass objects. 

Stellar mass black holes and active galactic nuclei 
(AGNs) appear to lie on a 'variability plane' defined 
by a relation between the black hole mass, accretion 
rate, and characteristic frequency in the PSD {yi, frc- 
quenc y of the lower high-fr e quency Lorent zian c omp o- 
nent) (|McHardv et al.l 120061: iKording et all 120071) . Due 
to the limite d sens itivity, v\ is not detected in ULXs. 
ICasella et al.l |2008) proposed a means to estimate the 
black hole mass in ULXs with type C QPOs via the vari- 
ability plane by assuming a relation between the QPO 
frequency and v\. Thus in the efficient accretion case, 
the black hole mass can be written as lg(M/M Q ) > 
0.51g(Lx/0.1c 2 ) — 0.51 lg^QPo — 7.9, where the inequality 
is due to the unknown bolometric correction. Given an 
X-ray luminosity of 10 40 erg s^ 1 and a QPO frequency 
of 8.5 mHz, the estimated mass is 1500 Mq or slightly 
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higher (due to bolometric correction). The 'variability 
plane' is tested for a wide range of black hole masses. 
However, the correlation between i/j and i'qpo na s not 
been tested in ULXs and the unknown bolometric cor- 
rection introduces additional uncertainty. 

The energy spectrum of the source from the two 
X MM-Newton observations is similar to that reported 
bv lRoberts fc Colbertl (|2003|) from Chandra data, best- 
fitted with a model consisting of a cool thermal disk com- 
ponent with kT rs 0.18 keV plus a power-law tail with 
a photon index of ~ 2.2. The best-fit inner radius of 
the accretion disk is several 10 3 km, which gives an up- 
per limit of ~10 3 Mq for a Schwarzschild black hole or 
~10 4 Mq for an extreme Kerr black hole. This is consis- 
tent with the value derived from the characteristic time 
scales. However, this mass estimate depends on inter- 
pretation of the soft component disk emission and that 
there is debate about the correct spectral model to apply 
to ULXs. 

To summarize, we have detected strong X-ray variabil- 
ity in NGC 6946 X-l, which becomes the most variable 
ULX known so far. A break in the broad band timing 
noise is firmly identified at a frequency near 3 mHz. A 
tentative signature of QPOs, likely of type C, is seen. 
The source is most likely in the intermediate state and 
has behavior analogous to that of NGC 5408 X-l and 
some states of GRS 1915+105. Several different rela- 
tions known to hold for stellar and/or supermassive black 
holes, each with their own caveats and systematic uncer- 
tainties, suggest that the characteristic time scales of the 
PSD imply the presence of a black hole with a mass on 
the order of 10 3 M Q in NGC 6946 X-l. 

We thank the anonymous referee for helpful comments. 
HF acknowledges funding support from the National 
Natural Science Foundation of China under grant No. 
10903004 and 10978001, the 973 Program of China under 
grant 2009CB824800, and the Foundation for the Author 
of National Excellent Doctoral Dissertation of China un- 
der grant 200935. 

REFERENCES 

Abolmasov, P., Fabrika, S., Sholukhova, O., & Kotani, T. 2008, 

arXiv: 0809.0409vl 
Begelman, M. C. 2006, ApJ, 643, 1065 
Bclloni T., Hasinger G. 1990, A&A, 230, 103 

Belloni, T., Homan, J., Casella, P., van der Klis, M., Nespoli, E., 
Lewin, W. H. G., & Miller, J. M., & Mendez, M. 2005, A&A, 
440, 207 

Belloni, T., Psaltis, D., & van der Klis, M. 2002, ApJ, 572, 392 
Blair, W., & Fesen, R. 1994, ApJ, 424, L103 



Blair, W. P., Fesen, R. A., & Schlegel, E. M. 2001, AJ, 121, 1497 
Casella, P., Belloni, T., & Stella, L. 2005, ApJ, 629, 403 
Casella, P., Ponti, G., Patruno, A., Belloni, T., Miniutti, G., & 

Zampieri, L. 2008, MNRAS, 387, 1707 
Colbert, E. J. M., & Mushotzky, R. F. 1999, ApJ, 519, 89 
Dewangan, G. C, Titarchuk, L., & Griffiths, R. E. 2006b, ApJ, 

637, L21 

Fabbiano, G., & Trinchieri, G. 1987, ApJ, 315, 46 

Feng, H., Kaaret, P. 2005, ApJ, 633, 1052 

Feng, H., Kaaret, P. 2007, ApJ, 668, 941 

Feng, H., & Kaaret, P. 2010, ApJ, 712, L169 

Feng, H., Rao, F., Kaaret, P. 2010, ApJ, 710, L137 

Fridriksson, J. K., Homan, J., Lewin, W. H. G., Kong, A. K. H., 

& Pooley, D. 2008, ApJS, 177, 465 
Harlaftis, E. T., & Greiner, J. 2004, A&A, 414, L13 
Herrero, A., Kudritzki, R. P., Gabler, R., Vilchez, J. M., & 

Gabler, A. 1995, A&A, 297, 556 
Heil, L. M., Vaughan, S., & Roberts, T. P. 2009, MNRAS, 397, 

1061 

Kaaret, P., Feng, H., Wong, D. S., & Tao, L. 2010, ApJ, 714, L167 
Karachentsev, I. D., Sharina, M. E., & Huchtmeier,W. K. 2000, 
A&A, 362, 544 

King, A. R., Davies, M. B., Ward, M. J., Fabbiano, G, & Elvis, 

M. 2001, ApJ, 552, L109 
King, A. R 2008, MNRAS, 385, L113 

Kording, E., Falcke, H., & Markoff, S. 2002, A&A, 382, L13 
Kording, E. G, Migliari, S., Fender, R. P., Belloni, T., Knigge, 

C, & McHardy, I. 2007, MNRAS, 380, 301 
Leahy, D. A., Darbro, W., Eisner, R. F., Weisskopf, M. C, Kahn, 

S., Sutherland, P. G, & Grindlay, J. E 1983, ApJ, 260, 160 
Markowitz, A., et al. 2003, ApJ, 593, 96 
Matonick, D. M., & Fesen, R. A. 1997, ApJS, 112, 49 
McClintock, J. E., & Remillard, R. A. 2006, in Compact Stellar 

X-ray Sources, ed. W. H. G. Lewin & M. van der Klis 

(Cambridge: Cambridge Univ. Press), 157 
McHardy, I. M., Koerding, E., Knigge, C, Uttley, P., & Fender, 

R. P. 2006, Nature, 444, 730 
Morgan, E. H., Remillard, R. A., & Greiner, J. 1997, ApJ, 482, 

993 

Mucciarclli, P., Casella, P., Belloni, T., Zampieri, L., & Ranalli, 

P. 2006, MNRAS, 365, 1123 
Muno, M. P., Morgan, E. H., & Remillard, R. A. 1999, ApJ, 527, 

321 

Papadakis I. E., & Lawrence A. 1993, MNRAS, 261, 612 
Poutanen, J., Lipunova, G, Fabrika, S., Butkevich, A. G., & 

Abolmasov, P. 2007, MNRAS, 377, 1187 
Roberts, T. P., & Colbert, E. J. M. 2003, MNRAS, 341, L49 
Schlegel, E. M. 1994, ApJ, 424, L99 

Soria, R., Motch C, Read, A. M., & Stevens, I. R. 2004, A&A, 
423, 955 

Shaposhnikov, N., & Titarchuk, L. 2009, ApJ, 699, 453 
Strohmayer, T. E., & Mushotzky, R. F. 2003, ApJ, 586, L61 
Strohmayer, T. E., & Mushotzky, R. F. 2009, ApJ, 703, 1386 
Strohmayer, T. E., Mushotzky, R. F., Winter, L., Soria, R., 

Uttley, P., & Cropper, M. 2007, ApJ, 660, 580 
van der Klis, M., 1989, in Timing Neutron Stars, ed. H. Ogelman 

& E. P. J. van den Heuvel (Dordrecht:Kluwer), 27 
Vignarca, F., Migliari, S., Belloni, T., Psaltis, D., & van derKlis, 

M. 2003, A&A, 397, 729 
Wijnands, R., & van der Klis, M. 1999, ApJ, 514, 939 
Zampieri, L., & Roberts, T. P. 2009, MNRAS, 400, 677 



